Global energy spectrum of the general oceanic circulation

Advent of satellite altimetry brought into focus the pervasiveness of mesoscale eddies \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{{{{{\bf{{{{{{{{\mathcal{O}}}}}}}}}}}}}}}}({100})$$\end{document}O(100) km in size, which are the ocean’s analogue of weather systems and are often regarded as the spectral peak of kinetic energy (KE). Yet, understanding of the ocean’s spatial scales has been derived mostly from Fourier analysis in small "representative” regions that cannot capture the vast dynamic range at planetary scales. Here, we use a coarse-graining method to analyze scales much larger than what had been possible before. Spectra spanning over three decades of length-scales reveal the Antarctic Circumpolar Current as the spectral peak of the global extra-tropical circulation, at ≈ 104 km, and a previously unobserved power-law scaling over scales larger than 103 km. A smaller spectral peak exists at ≈ 300 km associated with mesoscales, which, due to their wider spread in wavenumber space, account for more than 50% of resolved surface KE globally. Seasonal cycles of length-scales exhibit a characteristic lag-time of ≈ 40 days per octave of length-scales such that in both hemispheres, KE at 102 km peaks in spring while KE at 103 km peaks in late summer. These results provide a new window for understanding the multiscale oceanic circulation within Earth’s climate system, including the largest planetary scales.

Do the mesoscales of Oð10 2 Þ km actually form the peak of the ocean's KE spectrum? What is the KE content of scales larger than Oð10 3 Þ km which constitute the ocean's large-scale general circulation (e.g., gyres, western boundary currents, Antarctic Circumpolar Current (ACC)), and are directly coupled to the climate system? Below, we present the first KE spectrum over the entire range of scales resolved in data from satellites and 1/12°-resolution models at the ocean's surface, including the spectrum at planetary scales. We find that the spectral peak of the global extratropical ocean is at ≈ 10 × 10 3 km and is due to the ACC. We see vestiges of a similar peak in the northern hemisphere, which is arrested at a smaller amplitude and at smaller scales (≈ 4 × 10 3 km) due to continental boundaries. Another prominent spectral peak is at ≈3 × 10 2 km, and with an amplitude less than half that of the ACC. Yet, the cumulative energy in the mesoscales between 10 2 km and 5 × 10 2 km is very large (>50% of total resolved energy). We also report the first observation of a roughly k −1 power-law scaling over scales larger than 10 3 km in both hemispheres, consistent with a theoretical prediction from a quasi-geostrophic model forced by wind 13,14 , with the power-law scaling extending up to the ACC peak in the southern hemisphere.
Our results here open exciting avenues of inquiry into oceanic and climate dynamics, allowing us to seamlessly probe interactions between motions at scales Oð10 2 Þ km and smaller with planetary scales larger than Oð10 3 Þ km. We are able to do so using a coarse-graining approach developed recently to probe multi-scale dynamics on the sphere [15][16][17] .
Our methodology, described in the "Methods" section and in refs. 15,16, allows us to coarse-grain the ocean flow at any length-scale of choice and calculate the KE of the resulting coarse flow. By performing a 'scan' over an entire range of length-scales, we extract the so-called 'filtering spectrum' without needing to perform Fourier transforms 18 . The filtering spectrum and the traditional Fourier spectrum agree when the latter is possible to calculate, as demonstrated in ref. 18 and in the Methods section. Unlike traditional Fourier analysis within a box/subdomain, coarse-graining can be meaningfully applied on the entire spherical planet, including land/sea boundaries, and so allows us to probe everything from the smallest resolved scales up to planetary scales.
Given a velocity field u and a filter scale ℓ, coarse-graining produces a filtered velocity u ' that only contains spatial scales larger than ℓ, having had smaller scales removed (see Fig. 1 and the "Methods" section). Unlike standard approaches to low-pass filtering geophysical flows, such as by averaging adjacent grid-cells or block-averaging in latitude-longitude, the coarse-graining of ref. 16 used here relies on a generalized convolution operation that respects the underlying spherical topology of the planet, thus preserving the fundamental physical properties of the flow, such as its incompressibility, its geostrophic character, and the vorticity present at various scales. The KE (per unit mass, in m 2 /s 2 ) contained in scales larger than ℓ is While E ' quantifies the cumulative energy at all scales larger than ℓ, the wavenumber spectrum quantifies the spectral energy density at a specific scale, similar to the common Fourier spectrum. Following 18 , we extract the KE content at different length-scales by differentiating the coarse KE with respect to length-scale, where k ℓ = 1/ℓ is the 'filtering wavenumber' and { ⋅ } denotes a spatial average. Note that there is no factor of 2π in our definition of k ℓ . Ref. 18 identified the conditions on the coarse-graining kernel for Eðk ' ,tÞ to be meaningful in the sense that its scaling agrees with that of the traditional Fourier spectrum when Fourier analysis is possible, such as in periodic domains. The Methods section shows how the filtering spectrum agrees with the Fourier spectrum performed within an oceanic box region over length-scales smaller than the box. The filtering spectrum has the important advantage of quantifying larger  scales without being artificially limited by the box size and windowing functions, with such approaches used to synthetically periodize the data as required to use Fourier methods. Figure 1 visualizes the flow from both AVISO satellite data and NEMO reanalysis model data (see "Methods") from a single daily mean at scales larger than and smaller than 10 3 km, termed "gyrescale" and "mesoscale", respectively. The colour intensity illustrates the flow speed and is consistent with expectations that the largescale flow magnitude is primarily dominated by the western boundary currents, while scales smaller than 10 3 km are dominated by mesoscale fluctuations. In the upper panels of Fig. 1, we can see clearly several well-known oceanic gyre structures, including the Beaufort Gyre in the Arctic, the Weddell and the Ross gyres in the Southern Ocean near Antarctica, the subtropical and subpolar gyres in the Atlantic and Pacific basins, and the ACC. North Atlantic currents are also readily observable, including the North Atlantic Current, its northward fork to the Norwegian Atlantic Current, and the southward East Greenland Current. The agreement between AVISO and NEMO is remarkable.
It is worth emphasizing that the flows in Fig. 1 are derived deterministically from a single daily mean of surface geostrophic velocity data without further temporal or statistical averaging. Past approaches have used climatological multi-year averaging (e.g., refs. 19,20) or Empirical Orthogonal Function (EOF) analysis (e.g., refs. 21,22), which is a statistical approach that requires averaging long time-series. Coarse-graining allows us to derive the dynamics governing the evolution of the flow in Fig. 1 (e.g., ref. 15), which is not possible for EOF analysis, and to disentangle length-scales and time-scales independently and in a self-consistent manner to study interactions between different spatio-temporal scales that link large-scale forcing, the mesoscale eddy field, and the global-scale circulation.
Results and discussion Figure 2 shows the filtering spectrum for both the northern and southern hemispheres as obtained from Eq. (2) using surface geostrophic velocity data from both satellite altimetry and a highresolution model (see "Methods"). This is the first spectrum showing the oceanic energy distribution across such a wide range of scales, from planetary scales Oð10 4 Þ km down to Oð10Þ km.
The top and bottom panels in Fig. 2 plot the same spectrum in linlog and log-log scale, respectively. The top panel highlights the prominent spectral peak due to the ACC, which is more than twice the mesoscale peak. The bottom panel highlights the power-law scaling over different k ℓ bands. The inset in the bottom panel plots the same spectrum on a lin-lin scale and highlights the wide range of wavenumbers around 300 km.
Note the zero energy content at scales larger than Earth's circumference and that energy also decreases precipitously when approaching the smallest scales resolved by each of the datasets, both of which are physical expectations. It is not possible for simulation, satellite, or field data to capture all scales present in the natural ocean, which certainly has scales smaller than 100 km. There is excellent agreement between satellite data and the higher resolution model data used here down to scales ≈ 100 km, which indicates that all scales larger than 100 km are well-resolved by both datasets, whereas smaller scales (10-100 km) are reasonably resolved only in the model data.

Antarctic circumpolar current and oceanic gyres
Unlike previously reported KE wavenumber spectra using Fourier analysis on box regions (e.g., refs. 6,8,9), some of which show a peak at mesoscales O(10 2 ) km, our Fig. 2 reveals that the largest spectral peak occurs at scales approximately 100 times larger, at ≈ 10 4 km, and only in the southern hemisphere. Indeed, a circle of latitude at 50°S has a geodesic diameter of ≈ 8.9 × 10 3 km. This can also be seen from the yellow colour of the ACC in Fig. 1 (top panels), highlighting its contribution to KE at large scales. Additional support that this spectral peak is due to the ACC can be found in Methods, where we plot the zonally (east-west) averaged KE as a function of latitude at various scales larger than 10 3 km. We can see that the dominant contribution is from latitudes [60°S, 40°S], which are roughly the latitudes of the ACC. We also see a much weaker signal at latitudes [30°N, 40°N], which roughly aligns with the Gulf Stream and Kuroshio. Further corroborating our assertion, the spectral peak in the southern hemisphere seen in Fig. 2 has no analogous peak in the northern hemisphere. Figure 2 shows vestiges of a similar peak in the northern hemisphere, but this is arrested at a smaller amplitude and at smaller scales ( ≈ 4 × 10 3 km) due to continental boundaries.

Gyre-scale power law
The KE spectra from both hemispheres in Fig. 2 at scales larger than 10 3 km reveal a range of scales that exhibit a ∼ k −1 power-law. This scaling has been predicted by 13 (see also ref. 14) for baroclinic modes at scales larger than the barotropic deformation radius, but has not been observed until now. The barotropic deformation radius is about 2 × 10 3 km in the oceans 2 and the ocean flow tends to be surface intensified as expected in a baroclinic flow 23 . Thus, the k −1 scaling observed in Fig. 2 is consistent with 13 . Previous studies relying on Fourier analysis within box regions would have had difficulty detecting such scaling due to the box size artifacts. The k −1 extends to larger scales and peaks at scales ≈ 4 × 10 3 km in the north, which is broadly the scale at which the flow starts feeling continental boundaries and gyres form. This can also be seen from the bright yellow colour of the Gulf Stream and Kuroshio in Fig. 1, highlighting their contribution to KE at large scales. In the southern hemisphere, on the other hand, the k −1 scaling extends up to the scale of the ACC, which encounters no continental barriers (at the latitudes of the Drake Passage) as it flows eastward around Antarctica.

Mesoscale eddies
In Fig. 2, we find a second spectral peak between 100 and 500 km, centred at ℓ ≈ 300 km, that is associated with the mesoscale flow. While we can see from Fig. 2 that the mesoscales do not form the largest peak of the KE spectrum, their cumulative contribution between scales 100 and 500 km greatly exceeds that of scales larger than 10 3 km. This is because the mesoscale flow populates a wider range of wavenumbers compared to the gyre-scale flow (note the logarithmic x-axis and the inset on a linear x-axis). Indeed, integrating the energy spectrum in Fig. 2 within the band 100-500 km yields more than ≈ 50% of the total energy resolved by either satellites or the mesoscale eddying model in the extratropics. Coarse-graining allows us to determine this fraction of KE belonging to the mesoscales in the global ocean. This is because integrating the filtering spectrum over all k ℓ in Fig. 2 yields the total KE (as resolved by the data), which was not possible in past studies using Fourier analysis in regional boxes.
The power-law spectral scaling at mesoscales and smaller scales has been the focus of many previous studies (e.g., refs. 6,11,24,25) using Fourier analysis within box regions. While this is not our focus here, we observe that the overall mesoscale spectral scaling lies between k −5/3 and k −3 in Fig. 2, consistent with previous studies 24,26 . Note that mesoscale power-law scaling is more clearly seen in smaller regions (see Methods) as the mesoscale power-law and the corresponding wavenumber range change significantly depending on the geographical location (see Fig. 15 of ref. 11).

Characteristic velocity and energy content within key scale bands
From the spectra in Fig. 2, we partition the energy conservatively into four bands of interest: ℓ ≤ 100 km, 100-500 km, 500-10 3 km, and ℓ ≥ 10 3 km such that the sum of their energy equals total KE resolved in the data. From KE within a scale band, KE band , we can infer a characteristic root-mean-square (RMS) velocity, u rms = ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 2 × KE band p at those scales. The results are summarized in Table 1. The mesoscale band (100-500 km) has the highest RMS velocity, between 15 and 16 cm/s, and accounts for more than 50% of the total energy resolved in the data. Mesoscales are slightly more energetic (per unit area) in the NH than SH. There have been several past attempts to quantify the mesoscale fraction of total oceanic KE, which is often cited as being 80% 27 . However, such estimates have large uncertainties due to working in small representative box regions that do not cover the global ocean 28 or due to relying on subjective detection criteria for mesoscale eddies 3 or on a Reynolds (or temporal) decomposition of the flow 1 . Table 1 shows that all resolved scales smaller than 500 km constitute ≈ 90% of the total surface geostrophic KE in the global extratropical ocean. Table 1 allows us to also infer a characteristic advective timescale, τ meso = '=u rms = Oð25Þ days. The RMS velocity decreases significantly for larger scales, with hemisphere-asymmetries becoming more prominent. Within the ACC-containing band of ℓ > 10 3 km, the NH and SH RMS velocities are approximately 4.2 and 5.4 cm/s, respectively, and with an associated characteristic timescale τ gyre = '=u rms = OðfewÞ years.
Seasonality and spectral lag time Figure 3 shows the seasonality in surface KE as a function of lengthscale from both satellite and model data, which exhibit similar trends. The most striking feature of Fig. 3 is the approximately constant lag time between length-scales of the same ratio as they attain seasonal maxima (red) and minima (blue). Going from 10 km up to 10 3 km, length-scales that are ×2 larger experience a lag of ≈ 40 days in their seasonal cycle, such that in both hemispheres KE at 100 km peaks in spring while KE at 10 3 km peaks in late summer. Figure 3 plots the normalized deviation, or z-score, z(t) = (x(t) − μ)/σ, where μ and σ are the temporal mean and standard deviation of the spectrum x = EðtÞ at each scale k ℓ . A detailed regression analysis is in the Methods section. These results agree with and extend previous analysis within regional boxes 24,29,30 , which found that scales between 50 and 100 km have maximal KE in the spring while scales larger than 200 km (but smaller than the box) tend to peak with a delay of one to two months. Possible explanations for the seasonal variation in KE at different scales include the increased eddy-killing from winter's high winds 17 , and an inverse energy cascade from the submesocales which energizes mesoscales in spring months 24,30 . While Fig. 3 is suggestive of an inverse cascade, in which seasonal variations propagate up-scale at the rate we observe, it alone is not sufficient evidence (see ref. 31 for other possible causes) and a direct measurement of the cascade as in ref. 15 is required but is beyond our scope here.  Table 2 in "Methods" for uncertainty estimates.
Gyre-scales. At gyre-scales the surface flow is influenced directly by continental boundaries, wind, and buoyancy forcing. Indeed, at scales > 10 3 km in Fig. 3, there is a noticeable break in the seasonal trends we discussed in the previous paragraph. In the SH, where the ACC is not impeded by continents, we see from Fig. 3 a pronounced winter peak at ≈ 10 × 10 3 km, which correlates with maximal wind forcing 17 . Scales between 10 3 km and 3 × 10 3 km in both hemispheres peak in autumn, consistent with previous analysis showing an autumn maximum in the surface flow of western boundary currents 32-34 due to the upper ocean seasonal heating cycle. At NH scales larger than 5 × 10 3 km, the KE is too small to be meaningful (Fig. 2).

Outlook
Using both satellite and model ocean data, our coarse-graining characterization of surface geostrophic velocity revealed novel spectral features, including the ACC and basin-scale gyres. By partitioning kinetic energy across length-scales in a manner that conserves energy and covers the global ocean, we showed that length-scales ℓ < 500 km make an overwhelming contribution to surface geostrophic kinetic energy due to populating a wide range of wavenumbers, despite not forming the most prominent spectral peak. Based on prior characterizations of ocean energy 1 , we reason that these length scales are dominated by mesoscale features such as geostrophic turbulence, boundary currents, and fronts. Our analysis also revealed a characteristic lag time, with length-scales that are twice as large experiencing a lag of ≈ 40days in their seasonal cycle. The expanded spectral analysis spurs new questions and lines of inquiry. We hope future investigations will shed light on the dynamic coupling between the spectral peaks at the gyre scales and mesoscales, determine if the k −1 slope between the two peaks is indeed due to baroclinic modes 13,14 , and whether the characteristic spectral lag-time is caused by an inverse cascade.

Description of datasets
For the geostrophic ocean surface currents, we use a Level 4 (L4) postprocessed dataset of daily-averaged geostrophic velocity on a 1/4°grid spanning January 2010 to October 2018 (except for the seasonality analysis, where we use 2012-2016). The data is obtained from the AVISO + analysis of multi-mission satellite altimetry measurements for sea surface height (SSH) 35 . The product identifier of the AVISO dataset used in this work is "SEALEVEL_GLO_PHY_L4_MY_008_047" (https://doi. org/10.48670/moi-00148).
We also analyze 1-day averaged surface SSH-derived currents from the NEMO numerical modelling framework, which is coupled to the Met Office Unified Model atmosphere component, and the Los Alamos sea ice model (CICE). The NEMO dataset consists of weakly coupled ocean-atmosphere data assimilation and forecast system, which is used to provide 10 days of 3D global ocean forecasts on a 1/ 12°grid. We use daily-averaged data that spans four years, from 2015 to 2018. More details about the coupled data assimilation system used for the production of the NEMO dataset can be found in refs. 36,37. The specific product identifier of the NEMO dataset used here is "GLOBAL_MULTIYEAR_PHY_001_030" (https://doi.org/10.48670/ moi-00021).

Coarse-graining on the sphere
For a field ϕ(x), a "coarse-grained" or (low-pass) filtered field, which contains only length-scales larger than ℓ, is defined as where *, in the context of this work, is a convolution on the sphere as shown in ref. 16 and G ℓ (r) is a normalized kernel (or window function) so that ∫d 2 rG ℓ (r) = 1. Operation (3) may be interpreted as a local space average over a region of diameter ℓ centred at point x. The kernel G ℓ that we use here is essentially a graded top-hat kernel: We use geodesic distance, γ(x), between any location x = (λ, ϕ) on Earth's surface relative to location (λ 0 , ϕ 0 ) where coarse-graining is being performed, which we calculate using γðxÞ = R E arccos sinðϕÞ sinðϕ 0 Þ + cosðϕÞ cosðϕ 0 Þ cosðλ with R E = 6371 km for Earth's radius. In Eq. (4), A is a normalization factor, evaluated numerically, to ensure G ℓ area integrates to unity. A convolution with G ℓ in Eq. (4) is a spatial analogue to an ℓ-day running time-average. The above formalism holds for coarse-graining scalar fields. To coarse-grain a vector field on a sphere generally requires more care 16 , particularly for vector fields that need not be toroidal (2D non-divergent). However, as this work focuses on SSH-derived 2D non-divergent velocity fields, these concerns do not apply here. More details can be found in ref. 38.
Comparing coarse-graining to Fourier analysis. It is common to quantify the spectral distribution of ocean kinetic energy via Fourier transforms computed either along transects or within regions; e.g., refs. 10,12,39-42. This approach has rendered great insights into the length scales of oceanic motion and the cascade of energy through these scales [43][44][45][46][47] . However, it has notable limitations for the ocean where the spatial domain is generally not periodic, thus necessitating adjustments to the data (e.g., by tapering) before applying Fourier transforms. Methods to produce an artificially periodic dataset can introduce spurious gradients, length-scales, and flow features not present in the original data 18 . A related limitation concerns the chosen region size, with this size introducing an artificial upper length scale cutoff along with an artificial discreteness of wavenumbers that can bias the KE distribution across scales. In this manner, no scales are included that are larger than the region size even if larger structures exist in the ocean. Furthermore, the data is typically assumed to lie on a flat tangent plane to enable the use of Cartesian coordinates. However, if the region becomes large enough to sample the earth's curvature, then that puts into question the use of the familiar Cartesian Fourier analysis of sines and cosines. The use of spherical harmonics, common for the atmosphere, is not naturally suitable for the ocean, again since the ocean boundaries are irregular. These limitations mean that in practice, Fourier methods are only suited for open ocean regions away from boundaries, and over a rather limited regional size.
As a demonstration of both the validity and advantages of coarsegraining for energy spectra, consider Fig. 4. This figure reproduces the energy spectrum from Fig. 3 of ref. 6, and includes both the coarsegraining, and traditional Fourier energy spectra measured from the 1/ 12°NEMO dataset. Spectra are calculated for the 5°× 5°box centred at 164°E, 37°N, which corresponds to the Kuroshio extension.
For length scales ≲ 200 km, the three spectra generally agree very well, and all produce close to a −5/3 spectrum. Ref. 6 used a higher resolution dataset, and so the spectra disagree for scales < 20 km as expected. At larger scales, coarse-graining does not require tapering, and so the spectrum at scales ≳ 200 km is not contaminated by the shape of the tapering window or the box size. As a result, coarsegraining is able to detect that the spectrum for this region peaks at ≈ 250 km. An exact relation between Fourier wavenumbers k and filtering wavenumbers k ℓ is provided by eq. (16) in ref. 18, which shows that the filtering spectrum Eðk ' Þ at k ℓ is essentially a weighted average of the traditional Fourier spectrum E(k) over a range of Fourier wavenumbers k centred around k ℓ . Therefore, it is worth being mindful that agreement between filtering and Fourier spectra may not be perfect in practice, as discussed at length in 18 , due to the fact that compact spatial filtering kernels are not strictly local in k-space compared to a sharpspectral filter, which can lead to additional smoothing as a function of scale. This is the price paid for gaining spatially local information, such as allowing us to distinguish the tropics and the hemispheres (see also Figs. 1 and 5), since concurrently exact spatial and scale localization is forbidden by the uncertainty principle.
ACC as the spectral peak. In Fig. 5 we provide plots a visualization of the zonally-averaged kinetic energy for selected filtering scales. Scales larger than 10 3 km have a dominant contribution from latitudes [60°S, 40°S], roughly corresponding to the ACC, and another contribution over [30°N, 40°N], roughly corresponding to the NH western boundary currents. Scales larger than 5 × 10 3 km continue to show a clear ACC signal, with no NH signal since this filter scale is just beyond the NH gyre spectral peak. Finally, scales larger than 12 × 10 3 km have no distinct ACC signal, showing that the ACC has been fully removed by coarse-graining at this scale. Combined, these results provide further support for our claim that the 10 4 km spectral peak corresponds to the ACC.
Land treatment. When coarse-graining near land, it is necessary to have a methodology for incorporating land into the filtering kernel (c.f. ref. 38 for more in-depth discussion of land treatments). In the work presented here, we make the choice of treating land as zero velocity water. Since coarse-graining is essentially a 'blurring' (analogized with taking off ones glasses to have a blurrier picture), the land-water division itself also become less well-defined, and so treating land as zero-velocity water is both conceptually consistent and aligns with noflow boundary conditions. Additionally, this land treatment allows for a 'fixed' (or homogeneous) filtering kernel at all points in space, and as a result allows for commutativity with derivatives (e.g., divergence-free flows remain divergence-free after coarse-graining) 16 . Note, however, that only the true water area is used as the denominator when computing area averages (e.g., the NH area-averaged energy is the coarse KE summed over all NH cells, including land, divided by the water-area of NH). This is merely a normalization choice and does not affect our results. For NEMO data, the normalization area is 104 × 10 6 km 2 in NH and 155 × 10 6 km 2 in SH. For AVISO, the area varies with time due to sea ice coverage, between [90.5, 98.7] × 10 6 km 2 for NH and [141, 154] × 10 6 km 2 for SH.
Deforming kernel approach. An alternative choice is to deform the kernel around land, so that only water cells are included, at the cost of losing the homogeneity of the kernel. The benefit to this approach is that it does not require conceptually treating land as zero velocity water. However, it has the significant drawback that coarse-graining no longer commutes with differentiation and, as a result, does not necessarily preserve flow properties such as being divergence-free. Additionally, a kernel that is inhomogeneous (i.e., changes shape depending on geographic location) does not necessarily conserve domain averages, including the kinetic energy of the flow, and has the potential to both increase or decrease the domain average. More details are provided in ref. 38.
Comparing land treatments. Figure 6 presents the energy spectra, similar to Fig. 2 in the main text, using both deforming and fixed kernels for the single day 02 Jan 2015. The deforming-kernel spectra agree remarkably well with the non-deforming (fixed) kernel spectra, in that they present the mesoscales, ACC, and gyre peaks at similar scales. There are some quantitative differences, such as the deforming kernel SH spectra presents a slightly broader and higher-magnitude ACC peak.  However, at very large length scales information from the equatorial band and opposing hemisphere can become introduced through the filter kernel. To resolve this issue, we use a 'reflected hemispheres' approach, wherein one hemisphere is reflected and copied onto the other hemisphere, essentially producing a world with two north, or two south hemispheres. It is worth noting that reflected hemispheres and equatorial masking would not be necessary in a context where ageostrophic velocities are also considered and a global power spectrum is desired. They are used here because we wish to disentangle the power spectra of the geostrophic flow in each of the extra-tropical hemispheres. Figure S2 in the Supplementary Material (SM) shows the filtering spectra from NEMO without relying on hemisphere reflection, and is to be compared to Fig. 2 in the main text. The two are in qualitative agreement, with an ACC peak in the SH and mesoscale peaks in both hemispheres. Unsurprisingly, the spectra only deviate for very large filtering scales, where an increasing amount of extra-hemisphere information is captured by the large kernels. Specifically, the NH spectra has a third peak at scales ℓ > 10 4 km that is not present when using reflected hemispheres. This very large-scale peak is a result of the NH kernels capturing the ACC. It is worth noting, however, that the main ACC peak is still present in the SH spectra, as is the NH gyre peak at approximately ℓ = 3 × 10 3 km. Table 1 presented median values of the RMS velocity and percentage of total KE contained within various ℓ-bands. Table 2 presents the interquartile range (25th to 75th percentiles) to provide an estimate for the sensitivity of those values.

Seasonality
A 60-day running mean is applied to remove higher frequencies and allow us to better consider the longer-time trend, and individual years are averaged onto a 'typical' year for the purpose of comparison. Seasonality results use the 5 years spanning 2012-2016 for AVISO, and the 4 years spanning 2015-2018 for NEMO.
A useful statistical method for comparing signals is to compare the normalized deviations, or z-scores, of the signal. For a set of points where μ x and σ x are the mean and standard deviation of the x i . As a result, data points that are larger than the mean produce a positive z-score, while those smaller than the mean produce a negative z-score. Note that the normalized deviation (z-scores) in Fig. 3 are computed independently for each k ℓ , and so comparing magnitudes between scales is non-trivial.
Regression analysis of phase shift. Figure 3 presents a clear phase shift in the seasonal cycle as a function of length-scale. In order to quantify the phase shift, we need to first extract a meaningful set of (k ℓ , time) points. To that end, we extract, for each k ℓ , the (i) times corresponding to the lowest 10% (dark blue in Fig. 3), (ii) two middlemost 10% for the two zero-crossings (white in Fig. 3), and (iii) highest 10% (dark red in Fig. 3) of the normalized deviations presented in Fig. 3. This is essentially extracting the (k ℓ , time)-coordinates for the line of    Fig. S3. Vertical dashed lines separate the data sources, with the text along the top indicating the source data and hemisphere. Diamonds indicate the regression slope estimate, dark envelopes the 75% confidence interval, and light envelopes the 95% confidence interval. The r-value for each regression fit is printed alongside the corresponding slope distribution. The left-most illustration, separated by a thick black line, presents the mean (diamond) and confident intervals (envelopes) across the 16 regression analyses. In the legend, 'min' means the lowest 10% z-score, 'mid' the two middle-most 10% groups, and 'max' the highest 10%.
darkest red, darkest blue, and the two white lines, resulting in a total of four regression sets. Periodic phase adjustment to the days of the year is applied to maintain monotonicity in time, and the k ℓ grid is truncated to focus on regimes with a clear linear trend. The extracted data points are shown as the dots/vertical bars in Fig. S3 in the SM, along with their corresponding regression fits. Figure 7 presents the linear regression slope analysis for the data shown in Fig. S3. The different regression analyses generally agree well, with 12 of the 16 regression sets indicating a 35-45 day time-lag per octave of spatial scales. From this analysis, we conclude that lengthscales that differ by a factor of two (i.e., ℓ 1 /ℓ 2 = 2) have seasonal cycles that are off-set by 41 ± 3 days. Scales that differ a decade (ℓ 1 /ℓ 2 = 10) would correspondingly have a phase shift of 136 ± 10days, or roughly 4.5 months.

Data availability
All data needed to evaluate or reproduce the results in the paper are present either in the paper or is publicly available. Data underlying our analysis can be downloaded from CMEMS at https://marine. copernicus.eu/services-portfolio/access-to-products/.

Code availability
The coarse-graining analysis was performed using our code FlowSieve, which is publicly available at https://github.com/husseinaluie/ FlowSieve.